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Abstract 

(-») ; We propose a relativistic particle model without Grassmann variables 

q_) _ which, being canonically quantized, leads to the Dirac equation. Both V- 

<""j . matrices and the relativistic spin tensor are produced through the canon- 

ical quantization of the classical variables which parametrize the properly 
constructed relativistic spin surface. Although there is no mass-shell con- 
straint in our model, our particle's speed cannot exceed the speed of light. 
ff. ' The classical dynamics of the model is in correspondence with the dynam- 

ic. , ics of mean values of the corresponding operators in the Dirac theory. In 

£Nj " particular, the position variable experiences Zitterbewegung in noninter- 

f^**) ' acting theory. The classical equations for the spin tensor are the same as 

■ * . those of the Barut-Zanghi model of a spinning particle. 

o 

1 Introduction 

S^ ■ Starting from the classical works [1-4], a lot of effort has been spent 



on attempts to understand the behavior of a particle with spin on 
the basis of semiclassical mechanical models [5-16, 21-27]. One pos- 
sibility is to use the Grassmann (anticommutative) variables for 
parametrization of the spin space [11, 12]. In their pioneer work [11], 
Berezin and Marinov have suggested such a kind model. Their pre- 
scription can be briefly summarized as follows. For nonrelativistic 
spin, the noninteracting Lagrangian reads y(xj) 2 + §£i£i, where the 
spin inner space is constructed from vector-like Grassmann variables 
Ci; iiij = ~ £?'&■ Since the Lagrangian is linear on &, their conjugate 
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momenta coincide with £, 7Tj = |^ = i&. The relations represent 

the Dirac second-class constraints and are taken into account by the 
transition from the Poisson (Grassmann) bracket to the Dirac one, 
the latter reading 

{&>&}db = *<V (1) 

Dealing with the Dirac bracket, one can resolve the constraints, 
excluding the momenta from consideration. This gives a very eco- 
nomic scheme for the description of a spin: there are only three 
spin variables £j with the desired brackets (JTJ). Canonical quanti- 
zation is performed, replacing the variables by the spin operators 
proportional to the Pauli cr-matrices, Si = \o % ([<7\cr : >] + = 28^): 

[S i ,S j ] + =^-6 ij . (2) 

They act on the two-dimensional spinor space *$? a . Canonical quan- 
tization of the particle on an external electromagnetic background 
leads to the Pauli equation 

t 9$ / 1 /A e A , 9 . eh \ ^ 

' h m = W p - ~c A) - eA " - 2^ CTB ) *■ < 3) 

Relativistic spin is described in a similar way [11, 12]. 

The problem here is that the Grassmann classical mechanics rep- 
resents a rather formal mathematical construction. It leads to cer- 
tain difficulties [11, 16] in attempts to use it for describing the spin 
effects on the semiclassical level, before the quantization. Hence 
it would be interesting to describe spin on the basis of usual vari- 
ables. While the problem has a long history (see [5-13] and references 
therein), there appears to be no wholly satisfactory solution to date. 
It seems to be surprisingly difficult [13] to construct, in a systematic 
way, a consistent model that would lead to the Dirac equation in 
the the course of canonical quantization. It is the aim of this work 
to construct an example of a mechanical model for the one-particle 
sector of the Dirac equation. 

To describe the nonrelativistic spin using commuting variables, 
we need to realize the commutator algebra of the operators Si 

[Si,Sj]- = ihe ijk S k , (4) 



instead of the anticommutator one ([21) . This has been achieved in 
the recent work [15] starting from the Lagrangian 



S = dt ( —(ii) 2 H — Aiii + eA 

o„*2 1 
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+ Tr&i eijkUjBk) + — — + -r((wt) - a ) • (5) 



The configuration-space variables are #»(£), ojj(t), ^(t) and 0(t). 
Here Xj represents the spatial coordinates of the particle with the 
mass m and the charge e, Ui are the spin-space coordinates, g, (j> 
are the auxiliary variables, a = const and B = V x A. The second 
and third terms in Eq. (jSJ) represent minimal interaction with the 
vector potential Aq, A{ of an external electromagnetic field, while 
the fourth term contains the interaction of the spin with a mag- 
netic field. At the end, it produces the Pauli term in the quantum 
mechanical Hamiltonian. 

The Dirac constraints presented in the model imply [15] that spin 
lives on a two-dimensional surface of a six- dimensional spin phase 
space Ui, 7Tj. The surface can be parametrized with the angular- 
momentum coordinates Si = eijkWjftk, subject to S 2 = ^j-. They 
obey the classical brackets {Si, Sj} = e^^Sk- Hence we quantize 
them according to the rule Si — > Si. 

The action leads to a reasonable picture on both classical and 
quantum levels. The classical dynamics is governed by the La- 
grangian equations 

6 6 

mii = eEi + -e ijk XjB k S k diB k , (6) 

c mc 

g 
Si = — eij k SjB k . (7) 

mc 

We have defined E = -\$± + VA . Since S 2 w Ti 2 , the S-term 
disappears from Eq. ([6]) in the classical limit h — ¥ 0. Then Eq. fl6]) 
reproduces the classical motion on an external electromagnetic field. 
Notice also that in the absence of interaction, the spinning particle 
does not experience an undesirable Zitterbewegung. Equation (JTJ) 
describes the classical spin precession in an external magnetic field. 
On the other hand, canonical quantization of the model immediately 
produces the Pauli equation (J3J). 



Below, we generalize this scheme to the relativistic case, taking 
angular-momentum variables as the basic coordinates of the spin 
space. On this basis we construct the relativistically invariant clas- 
sical mechanics that produces the Dirac equation through the canon- 
ical quantization, and briefly discuss its classical dynamics. 

2 Relativistic spin surface 

Here we motivate our choice of the variables for describing the rel- 
ativistic spin. The dynamical model based on these variables is 
constructed in the next section. The relativistic equation for the 
spin precession can be obtained including the three-dimensional spin 
vector Si 01]) either into the Frenkel tensor $^, ^ v u v = 0, or 
into the Bargmann-Michel-Telegdi 4- vector S 11 , S^u^ = 0, where 
u^ is the 4-velocity vectoio Unfortunately, the semiclassical mod- 
els based on these schemes do not lead to a reasonable quantum 
theory, as they do not produce the Dirac equation through the 
canonical quantization. We now motivate as to how this can be 
achieved in the formulation that implies the inclusion of Si into the 
SO(2,3) angular-momentum tensor J AB of five- dimensional space 
A =(/x, 5) = (0,1,2, 3, 5), 77^ = (- + + + -). 

In the passage from nonrelativistic to relativistic spin, we replace 
the Pauli equation by the Dirac one: 

G^D" + mc)^(x^) = 0, (8) 

where p^ = —ihd^. The position of the particle is described in the 
standard way; the corresponding phase-space variables are x^, p u , 
{x»,p u } PB = rf v , Tf v = (- + ++). 

Let us look for the classical variables that could produce the 
T -matrices. According to the canonical quantization paradigm, the 
classical variables, say z a , corresponding to the Hermitian operators 
z a must be chosen to obey the quantization rule 

[z?,z f> ]-=ih{J*,z fi } . (9) 

z— >z 

In this equation, [ , ]_ is the commutator of the operators and { , } 
stands for the classical bracketo To avoid the operator-ordering 

lr The conditions $'"'«„ = and S^u^ = guarantee that in the rest frame, only three 
components of these quantities survive, which implies the right nonrelativistic limit. 
2 It is the Poisson (Dirac) bracket in a theory without (with) second-class constraints. 



problems, we will consider only the sets of operators which form 
the Lie algebra, [z a , z"}- = c°^ 7 i 7 . So our first task is to study the 
algebra of T -matrices. We note that commutators of T M do not form 
closed Lie algebra, but produce SO(l,3) Lorentz generators 

[T»,T u ]_ = -2iT» u , (10) 

where T^ = i(TT v -T v T lt ). The set V, V" forms a closed algebra. 
Besides the commutator (TTUj) . we have 

[r^, r°]_ = 2i(^ a r u - n va r^), 
[r^, r Q/3 ]_ = 2i(rf MX r vfi - ^v ua - n vo r^ + r^r^). (11) 

The algebra can be identified with SO (2, 3) Lorentz algebra with 
generators L AB : 

[L AB , L CD ]^ = 2i(r, AC L BD - V AD L BC - V BC L AD + V BD L AC ), (12) 

assuming that T» = i>, Y' w = U v . 

According to Eqs. ([H]) and (I12p we need classical variables with 
the algebra 

{ J AB , J CD } = 2( V AC J BD - V AD J BC - V BC J AD + V BD J AC ). (13) 

The problem is that in classical mechanics the basic phase-space 
variables, say u a , rrp, necessarily obey the Poisson bracket {u a , irp} = 
S a /3. The algebra differs from ([TBI . So we generally need to pass 
from the initial to some composed variables as well as to impose 
some constraints. This implies the use of the Dirac machinery for 
constrained theories [16, 17]. 

To arrive at the algebra (fT3"|) . we introduce, tentatively, the ten- 
dimensional "phase" space of the spin degrees of freedom, u A , n B , 
equipped with the Poisson bracket 

{co A ,7r B } PB = rj AB , (14) 

and define the inner angular momentum 

jAB _ r./, ,AR , ,BA\ /-, r\ 

J = Z(uj n — <jj % ). (15) 

The Poisson brackets of these quantities form the desired classical 
algebra (TT31 . 

Below we use the decompositions (i,j = 1,2,3) 

jAB = (jS^ j^ = ( J50 5 j5i = j5 jOi = W? jij = ^ijk jjk y ^ 



Since the J AB are the variables which we are interested in, we 
try to take them as coordinates of the space u A , % B . The Jacobian 
of the transformation (u A , ti b ) — > J AB has rank equal to 7. So, only 
seven among the ten functions J AB (u,n), A < B, are independent 
quantities. They can be separated as follows. By construction, they 
obey the identity 

e^j\j a p = 0, <S> J ij = (J 50 )- 1 ^ 5 ^ ' - J 5j J 0i ), (17) 

that is the 3- vector D can be presented through J 5 , W 

D = -^J 5 x W. (18) 

Hence we can take J 5fl and J 0t as a part of the new coordinate 
system. Suppose, we complete this set up to a basis of the phase 
space adding three more coordinates, say a, b, c. Quantizing the 
complete set we obtain, besides the desired operators L 5fl , L°\ some 
extra operators a, b, c. They are not present in the Dirac theory, 
and are not necessary for the description of the spin. So we need 
to reduce the dimension of our space from 10 to 7, imposing three 
constraints. There is one important restriction on the choice of 
constraints. Canonical quantization of a system with constraints 
implies replacement of the Poisson by the Dirac bracket; the latter 
is constructed with the help of the constraints. We need SO(2,3)- 
invariant constraints T a , {T a , J ab }pb = 0; otherwise the Dirac- 
bracket algebra will not coincide with those of the Poisson, (fl3|) . 

The only quadratic SO(2, 3) -invariants which can be constructed 
from u A , n B are co a ua, oj A 7Ta and ti a tia- So we restrict our model 
to living on the surface defined by the equations 

T 3 = n A n A + a 3 = (19) 

T 4 = u A u A + a 4 = 0, T 5 = u A n A = 0, (20) 

where 03, a 4 are some numbers. Our suggestion is to take the surface 
as the inner space for the description of the relativistic spin. 

The matrix — — ', ija-, 5 '^ has rank equal ten. So the quantities 

.f\ J°\ T 4 , T 5 , u\ (21) 



can be taken as coordinates of the space (oj a , 7r B ). The equation 
jab _ 2{uj A 'k b - u b tt a ) implies the identity 

J AB J AB = 8[(^) 2 (7T B ) 2 - (A,) 2 ] = 

8[(T 4 -a 4 )(T 3 -a 3 )-(T 5 ) 2 ], (22) 

then the constraint T 3 can be written in the coordinates (12TJ) as 
follows: 

r , = (■/")' + 8 <f +a ,. (23) 

8(i 4 — a 4 ) 

where J u are given by Eq. ( Tl6|) . Note that T 3 does not depend on 
u 5 . On the hyperplane T 4 = T 5 = it reduces to 

- 8a 4 T 3 = ( J AB f - 8a 3 a 4 = 0. (24) 

Eq. ([23]) states that the value of 50(2, 3) -Casimir operatoi| ( J AB ) 2 
is equal to 8a 3 a 4 . 

In the dynamical model constructed below, the equation T 3 = 
appears as the first-class constraint. It implies that we are dealing 
with a theory with local symmetry, with the constraint being the 
generator of the symmetry [20]. The coordinate u 5 is not inert 
under the symmetry, 5u 5 ~ {T 3 ,u 5 } =7^ 0. Hence u 5 is gauge 
non-invariant variable. 

Summing up, we have restricted dynamics of spin on the surface 
(|TU1) . (I2"0"j) . If (|2"Tj) are taken as coordinates of the phase space, 
the surface is the hyperplane T 4 = T 5 = with the coordinates 
J 5fJ, ,J°\u 5 subject to the condition (1241) . Since u b is gauge non- 
invariant coordinate, we can discard it. 

This implies that we can quantize J 5m , J° l instead of the initial 
variables u A , n B . Similarly to the case for the T -matrices, the 
brackets of the variables J 5At , J° l do not form a closed Lie algebra. 
The nonclosed brackets are 

{ J 5i , J^} = {J°\ J°i} = -2J* j , (25) 

where J ij = {J m )- l (J u J°i - J^ J 0i ); see Eq. flT@. Adding them to 
the initial variables, we obtain the set J AB = (J 5m , J° l , J 1 ^) which 
obeys the desired algebra (1T3|) . According to Eqs. (jUJ), (fT2l) . (ITBI) . 

3 In quantum theory, for the operators 1121 . Jilt we have: J AB Jab = 20fi . 



quantization is achieved by replacing J AB on the T-matriceaj. We 
assume that u A has the dimension of length; then J AB has the 
dimension of the Planck's constant. Hence the quantization rule is 

This implies, that the Dirac equation (JSJ) can be produced through 
the constraint 

T 2 = p^ J 5m + mch = 0. (27) 

Summing up, to describe the relativistic spin, we need a theory that 
implies the Dirac constraints (1191) . (120]) . (1271) in the Hamiltonian 
formulation. 



3 Lagrangian action and the canonical quanti- 
zation 

One possible dynamical realization of the construction presented 
above is given by the following d = 4 Poincare-invariant Lagrangian 



2e 3 
1 

a 



[x» + e 2 ^) 2 - (e 2 u 5 ) 2 



[x^ + e 2 uj^)uju - e 2 u 5 u 5 



crmch a a$ 



2u 5 2e 

.A 



-e 4 (u A u A + a 4 ), (28) 



written on the configuration space x M , u^ 1 , u 5 , e*, a, where e,, a are 
the auxiliary variables. The local symmetries of the theory are the 
reparametrizations as well as the following transformations with the 
parameter e(r) (below we have defined j3 = e 4 e + \e 4 e) 

8x? = 0, 5u A = /3u A , 5a = /3a, Se 3 = 0, 

Se 2 = -pe 3 + %3, <5e 4 = -2e 4 /3 - ( ^V. (29) 

The presence of local symmetries implies the appearance of the first- 
class constraints ( TT9|) . (|27|) in the Hamiltonian formalism. We point 
out that in the Berezin-Marinov model the Dirac equation is implied 
by the supersymmetric gauge transformations. So, the symmetry 
(T2"9~j) represents the bosonic analogue of these transformations. 

4 The matrices TP, r M " are Hermitian operators with respect to the scalar product 



Although there is no the constraint p 2 + m 2 c 2 = in our model, 
our particle's speed cannot exceed the speed of light. Indeed, equa- 
tions of motion for the auxiliary variables e 2 , e 3 read (xuj) = e 2 a^ 
x 2 +2e 2 {xu) + e 2 ,oj A u)A + (? 2 a3 = 0. They imply (x M ) 2 = — ef>a 4 — a 2 a 3 . 
Since we are dealing with a reparametrization-invariant theory, we 
assume that the functions x M (r) represent the parametric equations 
of the trajectory x l (t). Then the previous equation implies 

(fJHfiy^-^H (3o) 

if we take 03 > 0, a 4 > 0. 

Curiously enough, the Lagrangian ( 1281) can be rewritten in al- 
most five- dimensional form. Namely, after the change lj (x^^a, e 2 ) 



->■ (x^ 
reads 


, x 5 , e 2 ), where x^ = x M - 


- ^uA x 5 = — ^-ui 5 

a ' a 


, e 2 = 


-■e 2 + 


(?)', ^ 


L = 


2e 3 2e 3 (o; 5 )- 
e 4 (. 


, . 4x9 e-imcTi 
^ ] + & + 

■jJ A U) A + 04), 


e 3 a3(' 
2{x' 


j 5 ) 2 
; ) 2 


(31) 


where 


we have defined Dx A = 


x + e 2 oo . 









Canonical quantization. To confirm that the action ([28]) leads to the 
Dirac equation, we construct its Hamiltonian formulation. In the 
Hamiltonian formalism, the equations ( fl9l) . ( )20|) . ( l27l) appear as the 
Hamiltonian constraints. The constraint T 2 has vanishing Poisson 
brackets with all the constraints. The remaining constraints obey 
the Poisson-bracket algebra 

{T 3 , T 4 } = -4T 5 , {T 3 , T 5 } = -2T 3 + 2a 3 , 

{T 4 ,T 5 } = 2T 4 -2a 4 . (32) 

If we take the combination 

f 3 = T 3 + ^T 4 , (33) 

a 4 

the algebra acquires the form 

{T 3 , T 4 } = -4T 5 , {T 3 , T 5 } = -2T 3 + 2^T 4 , 

a 4 



5 The change is an example of conversion of the second-class constraints in the Lagrangian 
formulation [18]. 



{T 4 ,T 5 } = 2T 4 -2a4. (34) 

The only bracket which does not vanish on the constraint surface 
is {T^,^}. According the Dirac terminology [16, 17], we have the 
first-class constraints (|33|) . ( 127|) . and the second-class pair (1201) . The 
presence of the first-class constraints indicates that we are dealing 
with a theory invariant under a two-parameter group of local (gauge) 
symmetries, which has been discussed above. 

The auxiliary variables a, e 4 turn out to be subject to their own 
second-class constraints. Assuming that the constraints are taken 
into account by the transition from the Poisson to the Dirac bracket, 
the variables can be omitted from consideration [16, 17]. The Hamil- 
tonian in terms of the remaining variables reads 

H = ^-(^) (n\ A + a 3 ) + ^(p,J^ + mch) + 
2e 3 \ 7r° J 2tc° 

Ae2^e2 + A e3 7r e3 . (35) 

Here 7r ea , are conjugate momenta for e a and X ea are the Lagrangian 
multipliers for the primary constraints ir ea = 0. 

The constraints (TT9|) . (120]) can also be taken into account by using 
of the Dirac bracket. Since they represent SO(2, 3) -invariants, the 
Dirac brackets of the quantities J AB coincide with the Poisson one, 
Eq. ( iTBl . Hence we quantize the model according to Eq. ( 1261) . In 
the quantum theory, the first-class constraint fl27|) is imposed on the 
state vector. This gives the Dirac equation. In the result, canonical 
quantization of the model leads to the desired quantum picture. 

4 Solution to the classical equations of motion 

We now discuss some properties of the classical theory and confirm 
that they are in correspondence with those of the one-particle sector 
of the Dirac equation [19, 20]. 

Besides the constraints discussed above, the Hamiltonian ( 1HB]) 
implies the following equations (we use the notation (pu) = p^u^): 

e a = Ka, T^ea = 0, a = 2, 3; (36) 



^ = 2^ J5 "' P " = 0] (37) 
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-* / 5 \ ^ 5 

W = — I e -2tL ^ + ^- p ^ fc" = e 2 p" - 2e 4 w M ; 

e^ V 7T 5 I 7T 5 

UJ 5 = - (^f]\ 5 + ^{pcol e = ^)-2e 4 u\ (38) 
e 3 y n a J it* n b 

Here e^ = ^- (^ph) f 2 " The equations reflect the fact that we are 
dealing with a theory with local symmetries. Indeed, we note that 
these equations do not determine the Lagrangian multipliers A ea , 
which enter as arbitrary functions into solutions to the equations of 
motion for x M , u A , ir A . According to the general theory [16, 17], 
variables with ambiguous dynamics do not represent the observable 
quantities. For our case, all the variables except p^ turn out to be 
ambiguous. 

To construct the gauge-invariant variables with unambiguous dy- 
namics we first note that x M , p M , J 5fl and J^ u represent e -invariant 
quantities. Equations for the angular-momentum variables follow 
from (138]) 

j 5 » = — %J^p u , P v = %{rfJ hv -p v J h "). (39) 

7T 5 7T 5 

In three-dimensional notation, these equations read 
j 50 = -(Wp), j 5 = -p°W + D x p, 

W = p°J 5 - J 50 p- (40) 

The ambiguity remaining in Eqs. (j3"7|) . (1391) due to the factor ^f has 
a well-known interpretation, being related to the reparametrization 
invariance of the theory. We assume that, in accordance with this, 
the functions x M (r), J AB {r) represent the dynamical variables x l (t), 

J AB (t) in the parametric form. Using the identity —^p- = c^^r, we 
obtain the deterministic evolution for the gauge-invariant variables: 

J r i J t5m 

§- = c(j™yijs\ *L- = -2c(j 50 )-V^, 

at at 

^— = 2c(J 50 )- (p"^ - p v J^). (41) 

To find the trajectory x l (t), we take the equations ( 1371) . ( 139|) in 
the gauge ^ = 1 and note that they imply the following closed 
third-order equation for x^ 

•& -#& = m^pv. (42) 
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For any p 2 < its solution is given by 

777 Cu (l^ 1 I 

*"(r) = 4 - ^P^ + -£=* cosier + ^). (43) 

Since x°(r) in a reparametrization- invariant theory must be a mono- 
tonic function of r, we take the integration constant a = 0. Then 
the parametric equations (143]) imply 

x*(t) = 4 + — *+-^=sCOs(o;t + ^). (44) 

Po V-P 

The solution is a combination of the rectilinear motion and oscilla- 
tions with the frequency 

mhp° 

For the particular value p 2 = —m 2 c 2 , and when p % <C p°, it ap- 
proaches to the Compton frequency ^e 2 -- Hence the variable x l (t) 
experiences the Zitterbewegung in noninteracting theory. 

As we have discussed above, our particle's speed cannot exceed 
the speed of light for any value of p 2 . For the case p 2 > 0, the 
general solution to Eq. (j42p represents hyperbolic motion 



x"(t) = 4 T p^r + -= cosh( Jp 2 r + (ff). (46) 

2p yp 

The existence of such a self-accelerated solution for the Frenkel elec- 
tron has been recently observed; see [21]. 

The variables free of Zitterbewegung. Besides the centre of charge, 
x, in the Dirac theory we can construct the centre-of-mass (Pryce- 
Newton-Wigner) [6, 7] operator x in such a way that the conjugated 
momentum of x turns out to be the mechanical momentum for x. 
So the Dirac particle looks like a kind of composed system (this 
picture has been used by Schrodinger [1] to identify spin with inner 
angular momentum of the system). The classical analogue of the 
centre-of-mass operator in our model is the variable 

5 M = of + -^ J^p v . (47) 

2p z 
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It obeys the equation x (r) = — f^rrP^; then the centre of mass 

x l (t) moves along the straight line, ^- — %-. Note also that p^ 
represents the mechanical momentum of the x^ -particle. 

Let us take as the classical four-dimensional spin vector the Pauli- 
Lubanski vector 5^ = e^ val3 p v J a p. It has no precession in the free 
theory; 5^ = 0. In the centre-of-charge instantaneous rest frame, 

J 50 = const, J 5 = 0, (48) 

it reduces to 5° = 0, S = p x W. According to Eqs. (Hgj) . f lT8|) . 
only the part W of the angular-momentum tensor ffl6l) survives in 
the nonrelativistic limit. 

Comparison with the Barut-Zanghi (BZ) model. The BZ spinning 
particle [14] is widely used [22-27] for semiclassical analysis of spin 
effects. Starting from the even variables z a , where a = 1,2,3,4 is 
the 50(1,3) spinor index, Barut and Zanghi have constructed an 
even spin tensor S^ v = \iz^^ u z. We point out that in the gauge 
% = 1 our equations ( |37|) . ()39l) coincide with those of the BZ model, 
making the identifications J 5/i -H- i> M , J 11 " -H- S^ . Besides, our 
model implies the equations (^-s-) 2 = — 4i?, p^J 5 ^ + mch = 0. The 
first equation guarantees that the centre of charge cannot exceed the 
speed of light. The second equation implies the Dirac equation_|. 



5 Conclusion 

In this work we have constructed a semiclassical model (j2"8|) . ( 13 5 p 
which describes the one-particle sector of the Dirac equation. Al- 
though there is no the constraint p 2 + m 2 c 2 = in our model, our 
particle's speed cannot exceed the speed of light. Spinning degrees 
of freedom are described on the basis of a seven- dimensional surface 
embedded in the ten-dimensional phase space u A , -k a equipped with 
the Poisson bracket (fT4"l) . The surface is specified by the 50(2,3)- 
invariant equations ffT^l) . (|2"U1) . The angular-momentum variables 
J 5m , J° l can be taken as coordinates of the surface. Quantizing 
them in accordance with their Poisson-bracket algebra (TlBl) , we have 
produced both the r M -matrices and the relativistic spin tensor J^ u . 
The first-class constraint (|27|) is imposed on the state vectors, which 
leads to the Dirac equation. 



6 The Barut-Zanghi model [10] does not imply the Dirac equation. 
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Our model shows the same undesirable properties as those of the 
Dirac equation in the semiclassical limit. We have solved the clas- 
sical equations of motion and confirmed that the position variable 
x l experiences the Zitterbewegung in noninteracting theory, see (1441) . 
The variable x 1 * specified by the equation (1471) moves along a straight 
line and corresponds to the Pryce-Newton-Wigner operator of the 
Dirac theory. Like the Dirac equation, the model presented here 
gives no evidence as regards which of these two variables should be 
identified with the particle position. 

We finish with a brief comment on a modification which solves 
the problems. We recall that the Dirac equation (jSJ) implies the 
Klein-Gordon one. In contrast, in classical mechanics the corre- 
sponding constraint (127)) does not imply the mass-shell constraint 
p 2 + m 2 c 2 = 0. So, the model presented here is not yet in complete 
correspondence with the Dirac theory. The semiclassical model that 
produces both constraints has been discussed in the recent work [28]. 
The extra first-class constraint implies that we are dealing with a 
theory with one more local symmetry, with the constraint being 
a generator of the symmetry [29, 17]. This leads to a completely 
different picture of the classical dynamics. The variable x^ is not 
inert under the extra symmetry; 5x^ = /3p M , where /3(r) is the local 
parameter. Being gauge non-invariant, x M turns out to be an unob- 
servable quantity. The variable x^ of Eq. (j4"7j) is gauge invariant and 
should be identified with the position of the particle. Because p^ is 
a mechanical momentum for x M , the particle's speed cannot exceed 
the speed of light. In the absence of interaction it moves along the 
straight line. Hence the modified model is free of the undesirable 
Zitterbewegung. 
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